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Abstract. Infrared and optical spectra of SN 1998bu at an age of one year after explosion are presented. The 
data show evidence for the radioactive decay of ^''Co to ^^Fe, long assumed to be the powering source for the 
supernova light curve past maximum light. The spectra provide direct evidence for at least 0.4 solar masses of 
iron being present in the ejecta of the supernova. The fits to the data also show that the widths of the emission 
lines increase with time. Photometric measurements in the H-band show that the supernova is not fading during 
the observation period. This is consistent with theoretical expectations. 
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1. Introduction 

Type la supernovae are believed to be the most prodigious 
producers of iron group elements and given the explosive 
nature of the events also the best recyclers of such ma- 
terial. The burning of a large fraction of a white dwarf 
near the Chandrasekhar limit into ^^Ni provides the en- 
ergy for the explosion and subsequent radioactive decay 
of the nickel to cobalt and then to stable iron provides the 
energy for the light curve. Type la supernovae are popu- 
lar distance indicators as they exhibit, following a simple 
calibration procedure, remarkable uniformity (see Phillips 
et al. [TM5|l . However, our understanding of these objects 
remains limited f Leibundeut 1200 I|l . 

The observation of the radioactive decay and accurate 
direct measurements of the mass of ^^Ni manufactured in 
the explosion provides for a better underpinning of the 
theoretical scenario outlined above. Evidence for radioac- 
tive decay comes from a number of sources. The light curve 
shape foremost suggests an exponentially declining power 
source. Kuchner et al. ('1994') showed that the flux ratio 
of a group of [Felll] lines to [Colli] increased with time 
in accordance with the expectations of radioactive decay. 
Schonfelder et al. H1996|l report on what would be the 
most direct and convincing evidence for radioactive decay, 
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namely the direct detection of the 7-ray lines originat- 
ing from the decay. However, even for the very luminous 
SN 1991T the detections are only marginal, although con- 
sistent with the current theories. For SN 1998bu Georgii 
et al. H2002fl report upper limits for the amount of ^^Ni 
manufactured in the explosion from the absence of ^^Co 
7-ray emission. Their 2-a upper limit is 0.35 M© of ^^Ni. 

The determination of the mass of ^^Ni manufactured in 
the explosion based on optical and infrared spectroscopy 
has been the subject of many papers and almost as many 
models. In this paper we present data collected in three 
observing runs one year after the explosion of the super- 
nova 1998bu. Two epochs of infrared (IR) data are used 
to explore the evolution of the radioactive element abun- 
dances. The later epoch is combined with optical data to 
derive some limits on the mass of iron present in the ejecta 
of SN 1998bu. 

2. Observations 

Supernova 1998bu was discovered in NGC 3368 (Villi 
I1998f) on May 9, 1998 and reached maximum light in the 
i3-band on May 21st 1998. The supernova has been ob- 
served extensively by a number of groups (Suntzeff et al. 
HP^ Jha et ahUnSl Hernandez et al.'^nill. NGC 3368 
(M96) has a distance modulus of 30.25±0.19 measured us- 
ing Cepheid variables (Tanvir, Ferguson & Shanks ri999|l 
and a recession velocity of 900±50km s^^. The extinction 
towards the supernova is high with an Ay = 0.94 ± 0.15 
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Fig. 1. Spectrum of SN 1998bu obtained with SofI at the 
NTT on January 26, 1999. 

(Jha et al. 119991 see also Hernandez et al. l2()()()l who derive 
an Av of l.OiO.l). 

SN 1998bu was observed with the SofI infrared camera 
spectrograph at the 3.58-m ESO New Technology tele- 
scope at La Silla on January 26, 1999 at an age of 250 
days (past maximum blue light). Two grism settings were 
used providing a spectrum covering the J, H and K bands 
(Figure 1). The spectroscopic standard BS4903 was used 
to remove telluric features. The spectral resolution of the 
data is (A/(5A)=500. The absolute fluxing of the data was 
based on the photometry obtained on the same night with 
the same instrument in the H band. The photometric 
standards were P550C and S217D (Persson et al. [TMSjl . 
The supernova H-hsmd magnitude was 17.5±0.1. 

The ISAAC infrared camera spectrograph at the 8.2- 
m ESO UTl telescope of the VLT array at Paranal was 
used to observe SN 1998bu on May 1, 1999 at an age of 344 
days. Imaging observations were obtained with ISAAC to 
photometrically calibrate the spectrum. The photometric 
standard was S860D (Persson et al. 1998 ) and the derived 
iJ-band magnitude was 17.5±0.1. The spectrum obtained 
covers the 1.2-1.9-/im wavelength range (Figure 2). The 
spectral resolution of the data is {X/6X)—15Q0. The ob- 
servations were obtained in two wavelength settings and 
the spectra were joined without any scaling. The accuracy 
of merging spectra obtained in different settings without 
scaling can be confirmed during the fitting of the data. The 
1.644-/im and 1.257-/im lines of [Fell] originate from the 
same upper level and their ratio is fixed by the radiative 
transition rates and the extinction only. As will be shown 
below the fits are excellent adopting standard atomic data 
and the commonly assumed extinction for SN 1998bu. 

FORSl observations of SN 1998bu at the 8.2-m ESO 
UTl telescope of the VLT array were obtained on June 
6, 1999 at an age of 381 days (past maximum blue light). 
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Fig. 2. Spectrum of SN 1998bu obtained with ISAAC at 
UTl on May 1, 1999. 

An optical spectrum covering the 0.4-0.9 -/.tm wavelength 
range was obtained. The spectral resolution of the data 
is (A/(5A)=440. The spectrophotometric standard EG274 
was used to flux the spectra. The absolute flux of the 
spectrum was also checked against the photometry de- 
rived from the acquisition images. The night was non- 
photometric and we have used the local standards as es- 
tablished by Hernandez et al. f2000') (stars 2, 3, 4 and 
5) as well as star 6 of Jha et al. ( 1 999,) to derive the SN 
magnitude. The scatter in the 5 local standards was less 
than 0.05 magnitudes and no colour terms were applied. 
The F-band magnitude of the supernova was 19.96±0.10. 

3. Discussion 

3.1. The H-band light curve 

The observational prescription for IR light curves of type 
la supernovae is that they decline by 0.013 magnitudes 
per day and comes from the work of Elias & Frogel (4983|l 
and the more recent update by Meikle H2000|l . However, 
these works do not extend their coverage signiflcantly into 
the first year after explosion. Extrapolating this rule, su- 
pernova 1998bu when observed with ISAAC should have 
been 1.2 magnitudes fainter than observed with SofI. This 
is not consistent with our observations (see Figure 3). 

To double-check our zero-points and the accuracy of 
the photometry we performed photometry of 26 stars com- 
mon in the field of view for the two instruments. The cho- 
sen stars all had magnitudes within ±1.5 mag of the super- 
nova, and were not significantly contaminated by the host 
galaxy. For the stars of comparable brightness to the su- 
pernova we derive the same magnitude in both the ISAAC 
and SofI data with a scatter consistent with the photomet- 
ric errors in each observation. There is no systematic off- 
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Fig. 3. Light curve of SN 1998bu in the H-h&nd. Early 
data from Meikle H2000|l and template light curve from 
Elias (1983). For the observation at the epoch of 350 days 
the lower point (dotted circle) indicates the magnitude of 
the supernova after the removal of a possible contribution 
from continuum not arising from the ejecta. 



set between the two imaging data sets and it is clear that 
within the errors of the photometry the supernova has not 
faded between the two observations (±0.1 magnitude). 

The possibility that a fraction of the flux is coming 
from an echo (reflected light from earlier times that is 
contributing to the observed magnitude of the supernova) 
is examined below. However, the late emission line spec- 
trum is not consistent with a significant contribution from 
an echo. 

3.2. Removing the echo 

Cappellaro et al. H2001|l have demonstrated that the op- 
tical spectrum of SN 1998bu at an age of 600 days is 
dominated by an echo. HST observations have resolved 
a ring surrounding the supernova which is consistent with 
emission from an echo. Garnavich et al. H2()()l|l from HST 
data determine the magnitude of the echo to be 21.4 in 
the visible (a factor of 2 higher than that determined by 
Cappellaro et al.). We traced the echo spectrum from 
Cappellaro et al. and have smoothed and scaled it by 
the HST photometry. It is reproduced together with our 
FORSl spectrum in Fig. H At the epoch of the FORSl 
spectrum no significant continuum emission is expected 
from the supernova ejecta. The scaled spectrum fits well 
under the FORSl data and in all fits to the data we have 
removed this contribution from our data. 

The SofI data, obtained at a time when no contribu- 
tion from the echo was seen in the Cappellaro light curve, 
also show no evidence of a continuum. The ISAAC data, 
however, do show evidence for an underlying continuum 



at a level of 3xl0~^*erg s^^ cm^^A^^. Without a de- 
tailed model of the dust and its scattering properties it 
is difficult to establish the exact contribution to this con- 
tinuum by the echo. At the epochs under consideration 
in this paper, there is no evidence from other supernovae 
for a continuum due to emission from the ejecta. A level 
of 3xl0~^*erg s~^ cm~^A~^ is somewhat higher than a 
simple extrapolation of the echo spectrum scaled as men- 
tioned above. However, given the broader H-hand light 
curve it is reasonable to expect the IR echo to be some- 
what brighter than the optical one. We adopt this contin- 
uum level as a non-ejecta contribution. 

Adjusting the supernova magnitude for this continuum 
results in a H-band magnitude of 17.7 for the supernova 
ejecta at the epoch of the ISAAC observations. We adopt 
the continuum subtracted spectrum as the one originating 
from the ejecta. Our conclusions are only marginally af- 
fected if we have overestimated the contribution from the 
echo in the IR spectrum. 



3.3. Fitting the spectra 

As discussed in the introduction a number of authors have 
used the forbidden lines of iron in supernovae to derive 
masses in the ejecta. Due to a paucity of data at epochs 
of one year in the near infrared almost all studies have 
concentrated on the optical regime (see for example Ruiz- 
Lapuente & Lucy '1992j Ruiz-Lapuente et al. 119951 Liu, 
Jeffery & Schultz ll997| . Spyromilio et al. ifM?^ did anal- 
yse the spectra of SN 1991T covering a broad wavelength 
range. The optical regime is affected by the reddening to a 
significant extent and in the absence of IR data other au- 
thors have often also fitted the extinction. With the long 
baseline provided by the IR observations the constraints 
on the fitting are much stronger. However, the sensitivity 
of the longer wavelength transitions to temperature and 
electron density are lower. The combination of optical and 
IR data constrains the fit significantly. 

To combine the ISAAC and FORSl data into a single 
spectrum we have first removed the echo spectrum from 
the FORSl and ISAAC data following the prescriptions 
discussed above. We have then increased the flux in the 
FORSl data by 0.0138 magnitudes per day (measured on 
the published l/-band light curve, excluding the echo con- 
tribution) making up for the 41 days of difference in the 
epoch of the FORSl and ISAAC data. 

Our modeling is based on the techniques and mod- 
els described in Spyromilio et al. (,1992^) . The line identi- 
fications are the same as those of Liu et al. H1997|) for 
the optical lines and those of Spyromilio et al. I)1992|l 
for the infrared lines. The temperatures and densities for 
our models are of order 4000 K and 2xl0^cm^'^ for the 
SofI epoch data and the same temperature but lower elec- 
tron density (8xl0^cm~^) for the ISAAC epoch. For sim- 
plicity cobalt and nickel are treated in LTE. Since we 
only derive abundances from iron lines we treat these in 
NLTE. The atomic data are from Nussbaumer & Storey 
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Fig. 4. Spectrum of SN 1998bu obtained with FORSl at UTl on June 11, 1999. The lower spectrum is a trace of the 
late light echo spectrum (Cappellaro et al. i2001jl scaled by a factor of 2 upwards - see text for details 



H1988a|l . Nussbaumer & Storey lPM1|l . Garstang (fTH^ . 
Berrington et al. H1991|l . Grevesse, Nussbaumer & Swings 
ifWTji . Nussbaumer & Storey H1988b|l . 

From the clean 1.644 /xm line of [Fell] we can derive a 
mass of iron emitting given the emissivity from the NLTE 
calculations used above. Using a distance of llMpc for 
M96, Ay=0.94, and an observed flux from the 1.644 /im 
line of 1.7x 10"^'' erg s~^ cm~^ then the mass of singly 
ionised iron present in the ejecta is derived to be 0.4 Mq 
with a conservative error of 0.1 M©. The value depends 
solely on the emissivity for this transition. Doubly ionized 
iron has strong lines in the optical. However, the [Felll] 
lines are heavily blended and are of higher excitation mak- 
ing their emissivity less reliable. A good fit to the [Fe III] 
lines could be achieved with a variety of conditions. We do 
not believe that from our models a more accurate mass of 
Fe++ can be derived. We leave that to more sophisticated 
modeling attempts in a future paper. 

Contardo, Leibundgut & Vacca H2()0()|l have analyzed 
the light curves of a number of type la supernovae and find 
typically 0.5 Mq of ^^Ni to have been manufactured in the 
explosion. This is consistent with our result for SN 1998bu, 
but only marginally consistent with the claimed upper 
limit from absence of the 7-ray line emission (Georgii et 
al. 12002 , Milne et al. 12001. 

Fransson, Houck & Kozma 1)1996(1 predicted from the- 
oretical modeling that the near-infrared lightcurves would 
flatten out a year after explosion. As discussed earlier we 
observe this effect. The mass of iron we observe emitting is 
inconsistent with this being the classical infrared catastro- 
phe as described by Axelrod (,1980,1 whereby the flattening 



of the cooling curve combined with the exponential decline 
of the heating due to the radioactivity combine to pro- 
duce a rapid cooling of the ejecta resulting in only ground 
state energy levels being populated. Mid-infrared tran- 
sitions would remain observable but higher energy lines 
should vanish. In SN 1987A the infrared catastrophe was 
seen without the disappearance of the optical and near-IR 
lines. However, only a small fraction of the mass of iron 
was emitting in these lines. In our data presented here 
this is not the case. A significant amount of iron is still 
emitting in the optical and near-IR. 

The flattening of the light curve may be due to the 
gradual transition of the emission from the optical to the 
near-IR as the ejecta expand and the densities drop. The 
electron densities remain high enough that the ejecta con- 
tinue to be coupled to the heating source. 

We note that also the type la SN 2000cx displayed a 
constant late light curve in the near infrared. The flat 
lightcurve was observed in the J- and H-bands between 
360-480 days past maximum, and appears consistent with 
late time modeling (SoUerman, Kozma & Lindahl 120041 
SoUerman, Kozma, Lindahl et al. 12004)1 . 

3.4. Line widths and profiles 

From the fitting of the spectra we detect a change in 
the line widths of the emission lines. After correction for 
the recession velocity of the parent galaxy we find that a 
Gaussian profile fits the observed emission well on the red 
side of each line. The Gaussian profiles do not produce a 
good fit to the blue side of the emission lines (see Figs. 
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Fig. 5. A fit to the spectrum obtained with SOFI at the NTT using a simple NLTE [Fel], [Fell] and [Felll] model 
and LTE [Coll] The narrow (few pixel) spikes in the data are artifacts of the atmosphere and/or the detector. 



[7| and ISl . The supernova seems to emit more in the blue 
than in the red. 

Since the lines we are observing are forbidden and 
moreover, are emitted in a rapidly expanding medium, 
self absorption is not considered. However, a skewed dis- 
tribution of the emitting material could easily explain the 
line profiles. 

The velocities required for the fits were 
5000±200kms~i for the SofI data and 7000±200kms-i 
for the ISAAC data. The errors on the velocities are 
conservative as multiple lines are fit. This increase in the 
observed emitting volume implies that the outer regions 
of the ejecta increased their relative emissivity with 
respect to the centre between the two epochs. The excess 
emission in the blue part of the profile means that the 
width of the Gaussian is constrained by the red side of the 
line and the expansion velocity is a minimum. A broader 
skewed line profile would provide a better fit. Since all 
3 instruments were wavelength calibrated independently 
we exclude a calibration error as the source of the effect 
seen. 

The change of the expansion velocities observed sug- 
gests that the results of Mazzali et al. (|1998;1 relating the 



observed late-time expansion velocity to the energy of the 
explosion should be reconsidered as the measured velocity 
may depend on the exact epoch of observation. 

3.5. Direct evidence for radioactive decay 



Varani et al. H1990|l demonstrated in the case of SN 1987A 
that the near IR fines of [Coll] and [Fell] at 1.547 and 
1.533 /xm, respectively, could be used to determine the ra- 
tio of the mass of iron to cobalt emitting in the ejecta of a 
supernova. The ratio of the two lines changes as radioac- 
tive cobalt decays into stable iron. The two lines arise from 
the same ionization stage and their excitation potentials 
are very similar. The ratio of the two lines is very weakly 
dependent on temperature and it is expected that the ion- 
ization structure of the cobalt and iron will be very simi- 
lar. Therefore the ratio of the lines is expected to directly 
relate to the ratio of the masses of iron to cobalt. 

The following formula was derived by Varani et al. 
P99U|1 for the mass fraction: 
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Fig. 6. A fit to the combined FORSl and ISAAC spectrum using a simple NLTE [Fel], [Fell] and [Felll] model and 
LTE [Coll], [Coin], [Nill], [Mill]. The FORSl spectrum has been scaled to match the ISAAC one based on the 
photometric evolution (see text for details). The narrow (few pixel) spikes in the data are artifacts of the atmosphere 
and/or the detector. 



where T is the temperature of the gas in degrees Kelvin 
(4000 K in our case). From the fits to the data (see Figs. 
ElandlHl) we derive the ratios and plot them as a function 
of age in Fig. El 

Superimposed on tlic two data points is tlic change of 
the mass ratio of cobalt to iron as expected from tlie ra- 
dioactive decay. A pure ^^Ni initial mass is assumed for 
this plot with no contribution by ^^Ni. The optical cobalt 
and iron lines have been used in the past by other au- 
thors (see Kuchner et al. II994|I to derive similar results 
for type la supernovae. 



4. Conclusions 

We have used late time optical and infrared photometry 
and spectroscopy of the type la supernova I998bu to de- 
termine the physical conditions in the ejecta a year after 
the explosion. The mass of emitting iron is consistent with 
determinations by other methods and in agreement with 
theoretical predictions. We find that the supernova is fad- 
ing very slowly in the H-band one year after the explosion. 



We find that the emission lines originating within the 
ejecta appear broader and skewed at later times possibly 
indicating an assymetric distribution of iron within the 
ejecta. The widths of the emission lines also increases with 
time. The emission observed in the infrared [Fell] lines is 
consistent with 0.4 solar masses of Fe+ being present in the 
ejecta. The evolution of the spectrum in the near infrared 
shows direct evidence of the radioactive decay of ^^Co into 
s^Fe. 
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